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ABSTRACT 



Aims. We investigate the early evolution of two distinct populations of low-mass stars in globular clusters under the influence of 
primordial gas expulsion driven by supernovae to study if this process can increase the fraction of second generation stars at the level 
required by observations. 

Methods. We analyse N-body models that take into account the effect of primordial gas expulsion. We divide the stars into two 
populations which mimic the chemical and dynamical properties of stars in globular clusters so that second generation stars start with 
a more centrally concentrated distribution. 

Results. The main elfect of gas expulsion is to eject preferentially first generation stars while second generation stars remain bound 
to the cluster In the most favourable cases second generation stars can account for 60% of the bound stars we see today. We also find 
that at the end of the gas expulsion phase, the radial distribution of the two populations is still dilferent, so that long-term evolution 
will further increase the fraction of second generation stars. 

Conclusions. The large fraction of chemically anomalous stars is readily explainable as a second generation of stars formed out of the 
slow winds of rapidly rotating massive stars if globular clusters suffer explosive residual gas expulsion for a star formation efficiency 
of about 0.33. 

Key words, globular clusters: general - stellar dynamics - methods: N-body simulations 
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1. Introduction 

Globular clusters are self-gravitating aggregates of tens of thou- 
sands to millions of stars which have survived over a Hubble 
time. Many observations show that these objects are composed 
of (at least) two distinct stellar populations. The first evidence 
rests on the chemical analysis that reveals large star-to-star 
abundance variations in light elements in all individual clus- 
ters studied so f ar, while the iron a bundance stays constant 
(for a review see Gratton e t al ] l2004h . This includes the well- 
documented anticorrelations between C-N , O-Na, Mg-Al, Li- 
Na and F-Na jKrafti Il994t ICar retta et al. 2006, 2007, 2009; 



Lind et al 



Gratton et alJ 20071 : iPasquini et al...2007; .Bonifacio et al..,2007i: 



2009h . This global chemical patter n requires H- 

burning at high temperatures around 75 x 10^ K dArnould et alj 
Ll999; Prantzos et al. 2007). As the observed chemical pattern is 
present in low-mass stars both on the red giant branch (RGB) 
and at the turn-off which cannot have reached such high temper- 
atures, the abundance anomalies must have been inherited at the 
time of formation of these stars. 

Further indications for multiple populations in individuals 
GCs comes from deep photometric studies that have revealed 
multiple giant branches or main sequences. In o) Cen a blue 
main sequence has been discovered (Bedin et al. 200 4) which 
is presumably related to a high content in He dPiotto et alj|2005l: 
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IVillanova et al.ll20()7[). A triple main sequence has been discov- 
ered in NGC 2808 dPiotto et al.ll2007h . A broadening of the main 
sequence of NGC 6752 has also been discovered which can not 
be explained by binary stars (Milone et al. 2010). The additional 
blue sequences are explainable only by a higher He content of 
the corresponding stars which increases the opacity and shifts 
the effective temperature towards higher temperatures. He-rich 
stars are also the progenitors of blue horizontal branch stars 
seen in many globular clusters (see Caloi & D'Antona 20051 
l2007l) . Whereas no direct observational link between abundance 
anomalies and He-rich sequences has been found, theoretically 
this link is easily understood as abundance anomalies are the 
main result of H burning to He. 

These observed properties lead to the conclusion that glob- 
ular clusters born from giant gas clouds first form a genera- 
tion of stars with the same abundance pattern as field stars. 
Then a polluting source enriches the intracluster-medium with 
H-burning products out of which a chemically-different second 
stellar generation forms with a spread of chemical peculiarities. 
This scheme can explain at the same time the abundance anoma- 
lies in light elements and He-enrichment. 

Two main candidates that reach the right temperature 
for H-burning have been proposed to be at the origi n of 
the abundance anomalies (I Prantzos" & Charbonneil l2006h : (a) 
intermediate mass stars evolving through the thermal pulses 
along the asymptotic giant branch (hereafter TP-AGB), and 
(b) main sequence massive stars. After being first proposed by 
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ICottrell & Da Costal (Il98lh the AGB scenario has been exten- 
sively studied (Ventura et al."200r, '2002'; 'Ventura & D'Antona' 
2005a b e, 2008a b, 2009; Denissenkov & Herwig 2003; 
Karakas & Lattanzio 2003; Herwig 2004a bt iFenner et al.ll2004l: 
Bekki et al.,,2007; .Decressin et al.. .2009). In massive TP- AGB 



stars (M > 4 Mq), the abundance anomalies are supposed to be 
created at the bottom of the convective envelope through hot 
bottom burning. 

On the ot her hand, a s has been sugge sted by 

IWallerstein eTaP (Il987h and iBrown & Wallerstein' ('1993'), 
massive stars can also p ollute the inter-stellar medium (ISM) 
of a forming cluster (see [S mith 2006; Prantzos & Charbonnel 
12006) . In particular iDecressin et al. (2007b) show that fast 
rotating massive stars (with a mass higher than ~ 25 Mq) are 
good candidates for the self-enrichment of globular clusters. 
In the wind of fast rotating massive stars (WFRMS) scenario, 
rotationally-induced mixing transports H-burning products 
(and hence matter with correct abundance signatures) from 
the convective core to the stellar surface, and, provided initial 
rotation is high enough, the stars reach break-up velocity 
while on the main sequence. As a result, a mechanical wind is 
launched from the equator that g enerates a disk around the star 
similar to that of Be stars (e.g. iTownsend et alj l2004h . Later, 
when He-burning products are brought to the surface, the star 
has already lost a high fraction of its initial mass and angular 
momentum, so that it no longer rotates at the break-up velocity. 
Matter is then ejected through a classical fast isotropic radiative 
wind. From the matter lost through the disk, a second generation 
of stars may be created with chemical pattern in agreement with 
observations. 

In the present paper we mainly focus on the dynamical con- 
straints related to the WFRMS scenario. However both scenar- 
ios are facing a similar problem. Within standard assumptions 
(canonical IMF and conservation of the stars in the cluster), 
the amount of matter lost by the polluter stars is much smaller 
than the mass locked into the first generation low-mass stars. 
Indeed, if we assume that first generation massive stars follows 
the canonical IMlB stars in the mass range 25-120 Mq account 
for only 10% of the mass of the whole first stellar generation 
(see IDecressin et a l. 2007a) and their slow winds account for 
2.5% of the mass only. If pollution is due to AGB stars, a sim- 
ilar constraint arises: the wind released by stars between 5 and 
6.5 Mq, for which nucleosynthesis agrees with the observations 
according to Ventura & D'Antona (2008a), represents less than 
3%. After taking into account the possible dilution of these slow 
winds with the pristine gas present in th e ISM to explai n the 
observed Li abundance variation (Decre ssin et al.l l2007al) . we 
find that the mass available to form the second generation low- 
mass stars compared to the first generation of low-mass stars 
is only about 10%. This is in sharp contrast with observations, 
which show that more than half and up to 85% of the stars in 
GCs are second ge neration stars (Prantzos & Charbonnel 2006; 
ICarretta et al.l2008 ). i.e., which show anticorrelations in light el- 
ements. Thus a rather extreme reduction of the first generation 
stars relative to the second generation stars is needed to repro- 
duce the observations. However massive binaries have recently 
been p roposed as polluters of the proto-GC by de Mink et al. 
(120091) . In this case the mass-budget is more favourable as more 
slow winds are ejected and more second generation stars are 
formed. This could help to reduce the fraction between first and 



second generation stars found in the present paper therefore sup- 
porting a high-mass star pollution scenario. 

One possible way to reconcile the pollution scenario with 
the observations is to consider a top-heavy initial mass func- 
tion (IMF) of first generation stars. In the case of pollution by 
fast-rotating massive stars, an IMF slope as flat as 1.55 (com- 
pared to the canonical value of 2.3) is required to reproduce the 
high numbe r of stars with abunda nce anomalies in the cluster 
NGC 6752 (IDecressin et al.ll2007al) . wh ereas die AGB scenario 
requires an even flatter IMF slope (see iPrantzos & Charbonnell 
2006). 

A second way to reconcile the pollution scenario with obser- 
vations is to consider that first generation stars are preferentially 
lost from the cluster during its evolution so that an initially rel- 
atively small population of second generation stars can become 
the dominant population after several Gyr To allow this prefer- 
ential loss of first generation stars requires the GCs to be initially 
mass-segregated (i.e., that more massive stars occupy the central 
part of the clusters). In this case the matter released in the disks 
of massive stars is more concentrated in the cluster centre, and 
second generation stars are born in the centre while first genera- 
tion stars are present throughout the cluster. 

The viability of a self-enrichment scenario by fast-r otating 
massive stars has been recently explored by Decressi n et al.l 
( l2008l Paper I). They have shown that first generation low-mass 
stars are preferentially lost from the cluster, which is assumed to 
be initially in dynamical equilibrium and mass-segregated, be- 
fore two-body relaxation induces a spread of second generation 
stars and a full mixing of the cluster D'Ercole et al. (2008) find 
similar results with the AGB scenario. Afterwards, the evolution 
is smoother and the variation of the fraction of second generation 
stars takes longer. Any radial difference between first and sec- 
ond generation stars is erased after 10-12 Gyr of evolution as the 
cluster relaxation time (a few Gyr) is much shorter than the age 
of the clustersQ In Paper I we show that even if the relaxation- 
driven evaporation increases the fraction of second generation 
(which harbour abundance anomalies) to about 25%, this ratio 
rema ins too low to fully explain the observations (between 50- 
85%. ICarretta et aL200 9). The increase of the fraction of second 
generation stars mainly occurs in the early times and points to- 
wards the high sensitivity of the fraction of second generation 
stars on cluster dynamics. 

In this paper, we aim to quantify the increase of the frac- 
tion of second generation stars to the total number of low-mass 
stars by another dynamical mechanism not taken into account in 
the above studies, namely the effect of primordial gas expulsion 
(i.e., the fast ejection of the remaining gas left by star forma- 
tion after the onset of supernovae). Gas expulsion can strongly 
modify the total binding energy of the cluster and can lead to an 
efficient loss of first generation stars from the cluster We empha- 
sise that we discuss generic properties of gas expulsion models 
based on the simplified assumption that a cluster contains only 
two stellar generations with the same [Fe/H]. Multiple popula- 
tions with different [Fe/H] would require other physical mecha- 
nisms, whereby notably gas accretion form the surro unding inter 
stellar medium (Pflamm-Altenburg & Krouoa 2009) may play a 
role, and recycling of SN ejecta also ( Tenorio-Tagle et al. 20071). 
In § 2 we present the N-body models used in this study. Then 



' The canonical IMF is a two part power-law function, f(m) oz m 
with cfi = 1.3 for stellar mass es 0.08 < m/M f, < 0.5 and 02 = 2.3 
(Salpeter value) for m > 0.5 M^ jKroupall200lh . 



- The only exception is the GC a> Cen, for which the relaxation time 
at the center is comparable to its age. Indeed in this cluster stars on 
the blue main sequence (i.e.. He-ric h) are found more centrally concen- 
trated than red main sequence stars dVillanova et alj|2b07l : lBellini et al.l 
120091) . 
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our results are discussed in § 3. In § 4 we present a complete 
scenario for the evolution of GCs and our conclusions are in § 5. 

2. Description of analysis 

The results presented in t his paper are based on t he gri d of N- 
body models computed bv lBaumgardt & Kroupal (|2007|) . where 
the effects of primordial gas expulsion on the dynamics of star 
clusters were studie d. The N-body models are computed with 
the NBODY4 code (lAarsethlll99 9^ and follow the evolution of 
20000 single-mass stars. The gas is treated as a spherical addi- 
tional potential that is removed gradually in order to change the 
total binding energy of the cluster The initial cluster follows a 
Plummer distribution. The cluster evolution is computed for 100 
to 150 initial crossing times so that the cluster can settle into a 
new equilibrium configuration and two-body relaxation (which 
acts on a much longer timescale) is not an important parameter 
in these modeJS; 

iBaumgar dt & Kroupal ( |2007[) studied in particular the influ- 
ence of three physical parameters on the early cluster dynamics. 
The first is the star formation efficiency, e, given by the ratio 
between the stellar mass and the initial mass of the parent gas 
cloud. This parameter defines the fraction of gas converted into 
stars due to star formation. The second parameter is the ratio 
between the half-mass radius and the tidal radius, r^lr^, which 
quantifies the initial concentration of the cluster and the strength 
of the tidal field of the host Galaxy. Finally the third parameter is 
the ratio of the timescale for gas expulsion relative to the cross- 
ing time, tg£ /fci - This quantity determines the ability of stars to 
adjust their orbital parameters when the potential changes dur- 
ing gas expulsion. The full grid of models includes variation of 
the star formation efficiency between 0.05 and 0.75, of the ratio 
of the half-mass radius to the tidal radius between 0.01 and 0.2, 
and of the ratio of the timescale for gas expulsion relative to the 
crossing time betwee n and 10. 

As the models of 'Baum gardt & Kroupal (l2007h take into ac- 
co unt only one stellar population, we use the same method as 
in jPaper ll to split the stars into two populations according to 
their specific energy. The stars with the highest binding energy 
mimic second generation stars that are more centrally concen- 
trated, while the other stars are assumed to be members of the 
first generation. We choose an initial fraction of second gener- 
ation stars of 10% to be consistent with the pollution by fast 
rot ating massive stars in t he ca se of a canonical IMF slope. 
As iBaumgardt & Kroiipal (l2007l) used only single-mass stars 
we cannot study in detail the mass dependence of our results. 
However the gas expulsion process we investigate here (duration 
shorter than 10 crossing times) acts on a much shorter timescale 
than two-body relaxation (a few Gyr) which could lead to the 
preferentially loss of low-mass stars, and also acts on a shorter 
timescale than the lifetime of low-mass stars. So we do not 
expect our results to depend much on the mass of the stars. 
Similarly, the short duration of the gas expulsion phase com- 
pared to the two-body relaxation timescale allows us to infer 
results suitable for the study of the early dynamics of globular 
clusters even with the limited number of stars (20000) in the 
N-body model library. 

3. Results 

3. 1 . Analysis of individual models 

iBaumgardt & Kroupal (i2007i) find that gas expulsion can lead to 
all situations between a cluster totally disrupted and an unaf- 



fected cluster, depending on the parameter values for the gas- 
expulsion timescale, tqe, the star formation efficiency, e, and 
cluster concentration c. In this paper we will concentrate on in- 
termediate cases that predict a large stellar mass-loss but with a 
remnant dense core. In Fig. [T] we present three such interesting 
cases. 

Case 1. The left panels show predictions for the gas expul- 
sion parameters e - 0.2, c = rh/rt - 0.01, and TcE/ta - 1, 
which represent an initially concentrated cluster with a low star 
formation efficiency and a long timescale for gas expulsion. In 
the upper left panel we show the evolution of the number of stars 
of the first and second population still bound to the cluster We 
use the following criteria to define if, at a given time, a star is 
bound to the cluster: the star needs to be within the tidal radius 
of the cluster and to have a negative total energy (sum of the ki- 
netic and potential energy). Initially some stars have a positive 
total energy and they move away from the cluster However dur- 
ing about the first 10 crossing times the radius of the stars bound 
to the cluster expands so that these stars stay within the cluster 
tidal radius. During the first 30 crossing times, about 10% of the 
stars are in this situation (positive total energy and still within the 
tidal radius of the cluster) and it de pends on the criteria used if 
they are to be considered as bound (iBaumgardt & Kroupaii200'7l 
with only a radius criterium) or unbound (this study with both a 
radius and a energy criterium). However after around 30 cross- 
ing times, both criteria give the same results so that at the end 
of the computation no difference exists between both criteria. 
The two-phase decrease of the number of stars comes from the 
change of an energy dominating criterion in the early times to 
a radial criterion that dominates after 30 crossing times. Using 
only the radial criterion would ha ve lead to a smoother decrea se 
of the number of bound stars (see iBaumgardt & Kroupaii2007l) . 

As expected from the initial radial distribution, case 1 depicts 
a cluster that loses more of its first generation stars (about 80%) 
than second generation ones (about 12%) so that the fraction of 
second generation stars is around 30% at the end of the com- 
putation (middle left panel in Fig.[TJ. The loss of stars is very 
pronounced during the first few crossing times when the lower- 
ing of the cluster binding energy is driven by the gas expulsion. 

An interesting point is that the radial distributions (see left 
bottom panel in Fig.[T]) differ at the end of the simulation (about 
100 initial crossing times), second generation stars being still 
more concentrated than the first generation ones. Thus we can 
expect that the fraction of second generation stars will increase 
further due to the relaxatio n driven long-term evaporation of the 
clusters, as seen in'Paper 1 Moreover, the final radius of the clus- 
ter is much larger than initially. Indeed the half-mass radius of 
first and second generation stars is 1 .0 and 0.6 pc respectively at 
t -0, and becomes 8.5 and 3.8 pc at the end of the computation. 
This strong radial expansion by a factor 8.4 (first generation) and 
6.3 (second generation) is mainly due to the long timescale for 
gas expulsion, tce, that is similar to the crossing time, allowing 
stars to adopt new orbital parameters with wider orbits without 
being lost from the cluster It should be noted that the tidal field 
is also weak in this case as shown by the radial extension of the 
cluster up to about 30 times its initial half-mass radius. 

Case 2. The central panels of Fig. [T] present the case of a 
cluster near total disruption which loses ~ 95% of its first gen- 
eration stars due to the gas expulsion process. The initial param- 
eters are e = 0.30, r^/rt = 0.033 and TcE/ta- = 0.33. Compared 
to the previous case, this model has a smaller gas fraction af- 
ter star formation and gas expulsion occurs on a much shorter 
timescale. As this cluster shows a smaller increase of its radius 
during its evolution the two-phase decrease of the number of 
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Fig. 1. Top panels: number fraction of first (dashed lines) and second generation (full lines) stars relative to their initial number as a 
function of time. Each line is normalised to its initial number. Middle panel: fraction of second generation stars bound to the cluster 
as a function of time. Bottom panel: final (at 100 initial crossing times) radial distribution from the cluster centre for the stars of 
the first (dashed lines) and the second (full lines) generation. Right, central and left panels refer to three cases with different initial 
parameters indicated at the top. 



stars is limited to only the first 5 crossing-times. The number of 
bound stars goes through a minimum around 10 crossing times 
before increasing during the following 20 crossing-times. This 
behaviour is related to the strong ellipsoidal shape that the clus- 
ter displays during its expansion phase, leading to a significant 
number of stars lying in the outer part of the major axis of the 
cluster distribution where they are outside the tidal radius (and 
are hence counted as unbound stars). When the cluster contracts 
and becomes more spherical, part of these stars decrease their 



orbital radius below and become bound to the cluster again. 
At the end of the evolution, the cluster radius has only increased 
by a factor 3-4. The central part of the cluster is dominated by 
second generation stars (bottom panel) and the fraction of sec- 
ond generation stars at the end is above 60%. However due to the 
low number of bound stars in the simulation, statistics becomes 
too poor to precisely infer cluster properties. 

Case 3. Finally the right panels of Fig. [1] correspond to a 
model with initial parameters of e = 0.33, rh/rt = 0.06 and 
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Fi g. 2. Left: Fraction of s econd generation stars as a function of the final fraction of bound stars at the end of the computations 
of iBaumgardt & Kroupal ( l2007h . i.e., after about 100 initial crossing times. Dashed lines indicate limiting cases where no second 
generation stars are lost (upper) and no preferential loss of first generation stars occurs (lower). Estimates of the statistical errors are 
also included based on the number of first, A^i, and second, N2, generation stars bound to the cluster. Right: Final half-mass radius 
of the cluster as a function of the fraction of bound stars at the end of the computation. Clusters with different values of the star 
formation efficiency, e, are indicated with various symbols and colours. In both figures, numbered circles indicate the three models 
presented in Fig.[T]and discussed in § 3.1. 



Tcfi/fcr = 0.33 that also undergoes a strong loss of stars leading 
to a cluster with second generation stars counting for half of all 
cluster members. This case is particularly interesting as it has a 
small final half-mass radius which is only about twice the initial 
half-mass radius. Second generation stars dominate at the centre, 
thus we can expect that the further evolution of this cluster will 
increase the fraction of second generation stars to match the ob- 
served fraction (50-85%) of stars with anticorrelations in light 
elements when taking into account its whole evolution. 

3.2. The whole set of models 

Figure |2] (left panel) shows the fraction of second generation 
stars remaining bound after gas expulsion, /2G - N2l{N\ + N2), 
as a function of the fraction of remaining bound stars, /bound - 
(Ni +N2) /Mni ■ AH cases are located between two extreme scenar- 
ios: no loss of second generation stars (upper dashed line) and 
no preferential loss of first gener ation stars (horizo ntal dashed 
line). All the cases computed bv IBaumgardt & Kr oupa (2007) 
that have more than 100 bound stars at the end of the integra- 
tion are presented. To assess the statistical significance of the 

results we estimate the statistical error from cr = ^cr^ + cr^, 

where cri = A^j"''^ and 0-2 = are the statistical uncertainties 
for the number of first and second generation stars. In clusters 
which do not lose many stars (points on the right), the error is 
dominated by 0-2, while cri and 0-2 contribute significantly for 
clusters suffering a large loss of stars. For clusters with less than 
500 stars the remaining error is about 20%. 

A general trend is clearly visible in that the more dissolved 
clusters also have a higher fraction of second generation stars 



at the end. This behaviour is expected since second generation 
stars are initially more bound to the cluster and the gas expul- 
sion occurs on a timescale short compared to the two-body re- 
laxation timescale. For clusters near disruption, the fraction of 
second generation stars can be as high as 70-75% although with 
poor statistics. N-body simulations with a higher number of ini- 
tial stars will be very helpful to quantify these numbers more 
precisely. 

Even if the fraction of second generation stars is a monotonic 
function of the fraction of stars remaining bound to the cluster, 
other cluster properties present a more pronounced dependence 
with the initial parameters. Figure |2] (right panel) shows the fi- 
nal half-mass radius of the cluster as a function of the fraction of 
stars remaining bound to the cluster For clusters with high /bound 
the half-mass radius increases with decreasing value of the star 
formation efficiency. Indeed clusters with a high star formation 
efficiency are less perturbed when the primordial gas expulsion 
happens due to a lower mass of the gas remaining after star for- 
mation. This trend roughly remains at low e value but with more 
scatter. 

3.3. Initial conditions for globular cluster formation 

Figure [3] shows the fraction of second generation stars that re- 
main bound to the cluster after the end of the gas expulsion 
phase for all the input parameters used bv I Baumgardt & Kroupal 
(I2OOI . White areas indicate fully disrupted clusters while yel- 
low to red colours indicate the fraction of second generation stars 
at the end of the computation. By varying the SFE, we retrieve 
the three main behaviours for clusters. For clusters with low SFE 
(e between 0.05 and 0. 1), almost no cluster can survive. In con- 
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Fig. 3. Fraction of second generation stars, /2G, which remain bound after gas expulsion as a function of various initial parameters. 
Each panel has a different value for the star formation efficiency (from 0.05 to 0.75). White areas indicate dissolved clusters (i.e., 
with less than 100 bound stars). Black lines indicate levels where /2G is equal to 20% (outer contour) and 50% (inner contour). 
Blue dots in panels SFE=0.20, 0.30 and 0.33 indicate the three models presented in Fig.[T](Note that the ratio of the gas expulsion 
timescale to the crossing time is expressed in decimal logarithm). 



trast, gas expulsion has almost no effect for high SFE cases (see 
bottom right panel with SFE of e = 0.75). In this last case /2G 
remains always lower than 20%. 

Intermediate values of the SFE around 0.3-0.33 are more in- 
teresting as many cases lead to a high fraction of second gener- 
ation stars still bound to the cluster For a SFE of 0.33, a high 
fraction of second generation stars is obtained for concentrated 
clusters with a short timescale for gas expulsion (Tcslta <!)• 
These candidates could be good progenitors of real GCs (see 



second and third cases in Fig. [T]in § 3.2). It should be noted 
that a value of 0.33 for th e SFE is close to the one found by 
iParmentier & Fritz3 ( |2009|) from their study of the mass evolu- 
tion of c lusters and is also consistent with direct observational 
surveys dLada & L ad a 2003). As we have seen, clusters with a 
small rji/rt ratio are more prone to present a high fraction of sec- 
ond generation stars. Thus this fraction should be higher with 
large distance to the Galactic centre (larger rt). Alternatively, 
the tidal field may have been weaker because the Galaxy was 
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not ye t assembled. The observational trend observed bv lCarrettal 
(1200 6^ who shows that clusters with large orbital period and 
with high orbital inclinations relative to the Galactic plane pro- 
duce more extended O- Na and Mg-Al anticorrelations (see also 
iFraix-Burnet et al1l2009i) could be the imprint of the primordial 
gas expulsion process. 

In addition to a SFE around 0.33, other regions of parameter 
space favour a large increase of the fraction of second generation 
stars, namely a SFE around 0.25 combined with a fast timescale 
for gas expulsion and a very concentrated cluster (rh/rt < 0.05). 
A last possibility is for a higher SFE (e - 0.5) and an initially 
extended cluster (rh/rt > 0.1). However this case produces too 
extended clusters compared to the observed ones so that only 
SFEs around or below 0.33 are allowed to increase the fraction 
of second generation stars. 

4. Towards a complete scenario for the evolution of 
globular clusters 

As we have seen in the previous section, primordial gas expul- 
sion can be a very efficient mechanism to increase the ratio be- 
tween second and first generation stars. The most favourable 
physical conditions for globular cluster formation and early evo- 
lution are: (1) a star formation efficiency around e = 0.33, (2) 
a concentrated cluster relative to the tidal radius and (3) a fast 
timescale for gas expulsion. In the following we would like to 
consider how these constraints can be used to refine the scenario 
of the evolution of globular clusters with pollution by fast rotat- 
ing massive stars. 

4.1. The wind of fast-rotating massive stars scenario 

As already outlined in the introduction, our s cenario requires 
some b asic assumptions that are detailed in iDecressin et alJ 
(I2007ah . Let us recall here the main points. We suppose that the 
first stellar generation contains stars with initial masses between 
0.1 and 120 M© following a standard IMF with a Salpeter-like 
slope for stars more massive t han 0.8 Mp and a log-norma l dis- 
tribution for lower mass stars ( Paresce & De Marchill200 0'). but 
that second generation stars consist only of low-mass long-lived 
stars with initial masses between 0.1 and 0.8 M(£l We consider 
that the first generation polluters are fast-rotating massive stars 
(i.e., with initial masses above ~ 25 Mq) that enrich the ISM 
through their slow mechanical winds loaded with H-burning 
products. We assume mass segregation, primordial gas expulsion 
and long-term evaporation of first generation low-mass stars in 
order to reproduce the large fraction of second generation long- 
lived stars we see today (see § 3). 

The evolution of GCs passes through different key phases: 

1 . Formation of a first generation of stars. First generation stars 
(over the complete mass range 0.1 to 120 M©) form from a 
giant molecular cloud with a "normal" chemical composi- 
tion similar to that of contemporary halo field stars of similar 
metallicity. As shown in § 3, specific conditions are required 
at that phase: (1) an initial star formation efficiency for first 
generation stars defined as the total mass enclosed within 
first generation stars relative to the initial mass of the proto- 
cluster cloud in the same volume, eic, around 0.3 - 0.33 and 
(2) an initially highly concentrated cluster with a small half- 
mass radius (up to a few pc, see § 4.2.4). 

' The assumption about the mass range of second-generation stars is 
made only in order to minimise the constraints on the mass budget. 



2. Evolution of fast- rotating massive stars (m > 25 Mq) and 
cluster pollution. Part of the pristine gas that has not been 
consumed to form first generation stars must sit within the 
cluster during the lifetime of the less massive polluters, i.e., 
~ 7-10 Myr for a 25 Mq star (see § 4.2.4 below). Indeed the 
Li-free matter ejected by massive stars in their slow winds 
has to be mixed with pristine Li-rich intra-cluster gas in order 
to explain the Li-Na anticorrelation ob served in NGC 6 752 
(Pasquini et al. 2005J) a nd 47 Tuc (Boni facio et al.ll2007h and 
NGC 6397 (Lind et ani2009l) . This anticorrelation can actu- 
ally be used to constrain the amount of pristine gas involved 
in this dilution process (see § 4.2.1). 

3. Formation of second generation stars. Second generation 
low-mass stars (0.1 - 0.8 M©) form from the ISM mate- 
rial polluted to various degrees by the slow winds of mas- 
sive stars loaded with H-burning products. They have to be 
more centrally concentrated than their first generation coun- 
terparts, as requked by the number ratios b etween first and 
second generation objects we observe todav. IDecressin et al.l 
( 12007 a') propose that massive polluters of the first generation 
could be born in the centre of the cluster or could have mi- 
grated there rapidly through mass-segregation. In both cases, 
the second generation stars are created in their immediate 
vicinity and share a similar radial distribution. Here we de- 
fine the formation efficiency of the second stellar generation, 
€20, as the ratio of the total mass enclosed by second gen- 
eration stars to the total mass of the slow winds ejected by 
massive stars and the ISM matter used for the dilution pro- 
cess. 

4. Gas expulsion by SN. Then gas expulsion occurs and re- 
moves the interstellar gas left after the two episodes of star 
formatiorQ. For this process to be efficient enough the gas 
expulsion timescale must be very short, i.e., in the explosive 
regime, tge < fcr (see § 3). We propose that this process is 
induced by the supernova explosions of the first generation 
stars that did not contribute to the chemical pollution, i.e., 
with initial masses below 20-25 Mq (the more massive pro- 
genitors implode see § 4.2.3 below). During this phase most 
of the first generation stars that occupy the outer regions of 
the cluster are lost into the Galactic halo while second gener- 
ation stars are more centrally concentrated and remain bound 
to the GC. 

5. Long-term dynamical evolution. Later, the long-term 
relaxation-driven evaporation leads to the preferential loss 
of first generation stars over 2-3 relaxation timescales (see 
iPaper Ih . Finally, both populations are mostly mixed and 
no further evolution of the number ratio between first an d 
second generation long-lived stars is possible (Paper %. 
Nowadays only the globular cluster w Cen keeps a memory 
of the different initial distributions between first and second 
generations stars, because the two-body relaxation time in 
the core is comparable to the cluster age. 



4.2. Possible issues 

The scenario outUned above faces several issues that require fur- 
ther discussion. 



Note that in principle second generation massive stars could lead to 
the formation of a third generation of stars, but the higher order gener- 
ations would comprise only an insignificant fraction of the whole pop- 
ulation. 
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4.2.1 . Total star formation efficiency 

The first issue is related to the dilution process between ISM and 
the ejecta of the polluters, and more specifically to the amount of 
pristine gas consumed to form second generation stars. Indeed 
this process may change the total SFE, which is the main pa- 
rameter affecting the efficiency of gas expulsion as discussed in 

§3 3. 

iDecressin et a D (HoOTa) determined that in order to repro- 
duce the Li-Na anticorrelation in NGC 6752, the mass ratio be- 
tween pristine gas and slow stellar winds is around 1.15 after 
integration over time and IMF of the massive star polluters. On 
the other hand the mass lost by massive stars and recycled into 
the second generation represents only 3.5-4% of the total mass 
of the first generation stars. Thus for a star formation efficiency 
for the first generation of eic = 0.33, and assuming that all the 
matter ejected in the slow winds of the polluters is converted 
into stars with a dilution factor of 1 .2 with pristine gasQ we find 
that at most 1.32% of the protocluster gas is used to form the 
second sellar generation. Thus the SFE is only slightly modified 
(by a few percent) by the second episode of star formation. In 
other words, the SFE is determined mainly by the formation of 
the first stellar generation. 



4.2.2. Initial mass of proto-GC clouds 

We can now try to estimate the initial mass of the proto-gas 
clouds from which globular clusters are born. From the global 
value of the SFE (around 0.33) within the cluster forming vol- 
ume the proto-gas cloud is about three times more massive than 
the total mass of created stars of first and second generation. To 
explain the high fraction of anomalous stars obs erved in present- 
day glo bular clusters (aro und 85% in NGC 6752 Decressin et alj 
l2007aL 70% in NGC 2808 lPrantzos & Charbonnelll2006h a large 
fraction of the stars born in the cluster should have been lost 
from the cluster during its evolution. Decressin eLal- (2007a) es- 
timate that around 95% of first generation stars need to be lost 
in NGC 6752. Thus the initial mass of NGC 6752 should be at 
least 10 times more massiv e than its p resent-day mass. Given the 
luminosity of NGC 6752 (lHarry'1996) and a mass-to-light ratio 
of M/Lv - 3, we evaluate its actual mass to be ~ 3 x 10^ Mq. 
Thus the mass of the proto-gas cloud should have been of the 
order of 9 x 10*' Mq. 

However NGC 6752 is o ne of the most extre me cases in 
the number of anomalous stars. ICarretta et ai](l2008h statistically 
study 19 GCs and find that stars with abundance anomalies (their 
intermediate and extreme populations) represent 50 to 80% of 
cluster stars. Therefore NGC 6752 can be one of the most mas- 
sive GCs initially, and the initial mass of most GCs is only of or- 
der of several 10'' Mq. These high stellar masses would explain 
why the pre-superno ya feedback energy is n ot able to expel the 
gas from the cluster (iBaumgardt et alJl2008h . 




4 5 6 7 
log M^,^, (Mq) 
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Fig. 4. Top: binding energy of a gas cloud as a function of its 
mass with an initial half-mass radius of 0.5 pc (full line) and 3 pc 
(long-dashed lines). The total energy released by SN for stellar 
progenitor masses between 8 and 25 Mq is indicated as the short- 
dashed line. Middle: number of SNe needed to unbind a cluster 
as a function of its total mass with an initial half-mass radius of 
0.5 pc (full line) and 3 pc and the total number of SNe created 
by the cluster (in the mass range 8-25 Mq). Bottom: Crossing- 
time of a cluster as a function of its total mass with an initial 
half-mass radius of 0.5 pc (full line) and 3 pc (long-dashed line) 
and correspondin g gas expulsion timesc ale (full line and short- 
dashed line, e.g.. iBaumgardt et"ani2008l) . 



4.2.3. Detailed chronology: early cluster evolution 

One key issue concerns the necessity to retain pristine gas within 
the cluster during the lifetime of the polluters in order to account 
for the presence of lithium in the atmosphere of second genera- 
tion stars. This constraint is relevant for the issue whether very 
massive stars which produce the first SN are unable to clear out 
the cluster from its gas. 



' These assumptions are needed to maximise the initial number of 
second generation stars. 



This condition can also be fulfilled if at the end of their 
lives the polluter stars do not undergo supernova explosions 
but rather directly collapse into black holes, avoiding injec- 
tion of SN kinetic energy into the ISM. Observational clues 
from black hole X-ray binari es (iPortegies Zwart et aP 1 1 9971: 
lErgma & van den HeuvelllI998h as well as nucleosynth esis con- 
straints (lMaededlI992l: iKobulnickv & Skillmanlll997h suggest 
that black holes form from stars with masses above ~ 25 Mq 
(see also Hegeretal. 2003). This rough limit is confirmed by 
2D core-collapse SN simulations (lFrveiilI999h . which also pre- 
dict that progenitor stars more massive than ~ 40 Mq are unable 
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Fig. 5. Top: supernova rate as a function of time, using the 
iHurlev et al.k 2000) stellar-evolution routines, for a single stellar 
population of 10^ M© following a canonical IMF. Full and dot- 
ted lines refer to cases where all stars above 8 Mq produce a SN, 
or only stars in the mass range 25-8 Mq, respectively. Bottom: 
integration over time of the number of SNe. In both panels time 
is shifted to the SNe produced by stars with an initial mass of 

25 Mq. 



to launch shocks and thus do not produce a SN explosion. In 
view of the sensitivity of core collapse simulations to the (un- 
certain) input physics, these numbers must be considered with 
caution. For example, lowering the mean neutrino energy by 20 
% lowers the fallback black hole limit to ~ 15 Mq ( Frveiil 19991) . 
In addition, rotation is expected to decrease the lo wer mass 
limit for black hole formation (e.g.. lHirschi et al.l2004l) . Besides, 
iGeorgy et al.l (12009 ) found that the lowest masses which allow 
black hole formation decrease from 40 to 30 Mq when the metal- 
licity decreases from solar to Z = 0.004. 

Despite the theoretical uncertainties regarding the formation 
of black holes and the physics of SN explosions, what matters at 
this level of the discussion is the fact that the lower mass limit for 
a star to collapse directly into a black hole is very close to that of 
the fast-rotating massive polluter stars. This provides a natural 
way to avoid the deposition of SN kinetic energy into the ISM 



for the first 7-10 Myr of GC evolution, and thus to retain pristine 
gas within the cluster during the lifetime of the polluters. 



4.2.4. Detailed chronology: onset of gas expulsion 

On the other hand, if all massive stars with initial mass below 
25 Mq end up as SN we still have to check if they are numer- 
ous enough to expel the remaining pristine gas from the clus- 
ter potential well. Fig. |4] (top panel) compares the binding en- 
ergy for clusters with initial half-mass radii of 0.5 and 3 pc as- 
suming a SFE of 0.33 with the energy released by all SNe for 
stars between 25 and 8 Mq as a function of the cluster initial 
mass (gas and stars). T hese quantities are computed following 
Baumgar dt et alj (l2008l) . In the middle panel, we show the num- 
ber of SNe needed to unbind the cluster as well as the total num- 
ber of SNe produced from stars in the mass range 25-8 Mq. 

Both the binding and SN energies increase with the cluster 
mass but with different rates: the binding energy scales as Mj^^ 
while the SN energy scales linearly. For clusters with mass be- 
low 10^^ Mq (i.e., the mass of proto-GC gas clouds), SN from 
stars in the mass range 8-25 Mq produce enough energy for gas 
expulsion to operate. 



4.2.5. Fast gas expulsion 

Another issue concerns the timescale for gas expulsion, which 
has to be fast compared to the crossing time according to our 
scenario. 

To determine the timescale of gas expulsion by massive stars 
in the range 8-25 Mq we present in Fig.|5]the SN rates we expect 
from a stellar population following a canonical IMF (Kroupa| 
200 ih normalised to a 10'' Mq star cluster. We also indicate the 
total number of SNe the cluster has over time (bottom panel). 
After 30 Myr around 5000 SNe have exploded and even within 
1 Myr after the 25 Mq stars explode about 1000 SNe are pro- 
duced. 

We can now compare the time at which enough SNe are pro- 
duced to unbind a cluster to the crossing time of the cluster. This 
result is shown in Fig. |4] (bottom panel) where both timescales 
are indicated. For clusters with small half-mass radii (around 
0.5 pc), only clusters with initial mass below a few 10^ Mq ex- 
perience fast gas expulsion as needed. However such a mass is 
too low to be the initial mass of GCs. On the other hand, clusters 
with large half-mass radii (around 3 pc) can produce fast gas ex- 
pulsion up to a mass of 5 x 10^ Mq, which is consistent with the 
value we infer for most proto-globular clusters. In the extreme 
case of NGC 6752 (with mass up to 9 x lO*" Mq, see § 4.2.2) the 
initial half-mass radius should be of the order of 3-5 pc to allow 
a fast enough gas expulsion. 

This confirms the results of lParmentier & Fritzd (l2009h . who 
show that for clusters with mass around 10* - 10^ Mq gas ex- 
pulsion is more likely to happen in 2 or 3 crossing times im- 
plying an adiabatic regime (i.e., small radius ). B esides our find- 
ings are m ainly compatible with the ones of iMarks et al.l (120081 
Isubmited l) which explain the relation between the slope of the ac- 
tual mass function and the concentration of globular clusters as 
the dynamical response for the gas expulsion process. However 
the initial mass and radius of the cluster found by Marks et al] 
( Isubmitedl) are smaller than the ones obtain in this paper as they 
deduced a top-heavy IMF in order to obtain gas expulsion. 
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4.2.6. Place of birth of second generation stars 

Let us now address the issue of how the slow winds of fast ro- 
tating massive stars are recycled into second generation stars. 
[Decressin et al. (2007b) assume that the matter inside the equa- 
torial disc can already start to condense and produce a proto- 
stellar object. However this local formation of second gener- 
ation stars cannot allow the formation of a d istinct main se- 
quence as observed in u) Cen (Bedin et al. 2004) and NGC 2808 
(Pjot to et al. 2007) as pointed out by Decressin et al. (2007a) and 
iRenzinil (l2008h . 

However it is possible that the strong radiation pressure ac- 
celerates the disc so that it dissipates on larger scales. Here we 
assumed that the disc is dense enough so that stellar radiation is 
not able to accelerate the matter above the escape velocity of the 
cluster. The matter originating from the disc will be stored out- 
side the cluster centre and will fall back when it is cold enough. 
The main difficulty of this scenario is to prevent the mixing be- 
tween SN ejecta with these slow winds. This can be the case if 
the ISM (and the diluted slow winds) present inhomogeneities 
such that SN ejecta escape the cluste r mainly through low- 
density regions creating tunnels (see Prantzos & Charbonnel' 
[2006). A similar view is presented by Palous et al. (2009): ther- 
mal instabilities developed when the energy deposition to the 
ISM is dominated by stellar winds (with H-burning products). A 
large part of the matter sinks into the cluster centre in the form 
of compact high-density and cold gas which can be used to form 
second generation stars. Latter when SNe dominate the input en- 
ergy, the cluster is cleared out by a stationary outflow (see also 
iTenorio-Tagle et al.ll2007HWunsch et alj|2008h . 

In this case it will be possible that the slow component 
(i.e. the slow winds enriched with H-burning matter) remains 
decoupled from the SN ejecta and cools down so that it migrates 
towards the cluster centre. There the pressure of the ISM in- 
creases and second generation stars form. In this case second 
generation stars will have a much more homogeneous chemi- 
cal composition compared to the stochastic formation process 
near their massive progenitors. Therefore this scenario may re- 
produce the discrete He-sequences inferred from observations. 
However the continuous O-Na distribution found in many clus- 
ters is still challenging (see e.g.. .Carretta et al...2009.) . 

5. Conclusions 

In this paper we have studied the influence of primordial gas ex- 
pulsion by supernovae during the early dynamical evolution of 
globular clusters, in the context of clusters with two chemically 
and dynamically distinct stellar populations. In particular we in- 
vestigate if this dynamical process can explain the high number 
of observed second generation stars which harbour abundance 
anomalies in light elements. We deduce the following: 

- If the two populations have a diff'erent radial extent with 
second generation stars more concentrated, primordial gas 
expulsion is able to expel most of the first generation stars 
while most second generation stars can be retained. 

- For a given fractional mass loss by the cluster, the fraction 
of second generation stars is nearly independent of the gas 
expulsion parameters (see § 3.2). 

- The final observed fraction of second generation stars can 
constrain the initial properties of GCs as this fraction is high- 
est for clusters with SFE around 0.3-0.33, with concentrated 
clusters relative to the tidal field, and with a fast timescale 
for gas expulsion relative to the crossing time. 



- We infer proto-GC cloud masses of several 10* and up 
to 9 X 10* Mq for clusters which show a large fraction of 
chemically different second generation stars like NGC 6752. 
Their initial half-mass radii are in the range of ~ 1-3 pc (4- 
5 pc for the most massive cases). 

- It is possible to reproduce the fraction of second genera- 
tion stars in present-day GCs through cluster dynamical pro- 
cesses by combining gas expulsion and tidal stripping during 
long-term evolution of initially mass-segregated clusters. 

- The primordial gas expulsion process can a lso be at the ori - 
gin of the observational trend observed by Carrettal (l2006l) . 
who shows that clusters with large orbital period and with 
high orbital inclinations relative to the Galactic plane pro- 
duce more extended O-Na and Mg-Al anticorrelations. 

Acknowledgements. T.D. and C.C. acknowledge financial support from the 
french Programme National de Physique Stellaire (PNPS) of CNRS/INSU, and 
from the Swiss National Science Foundation (FNS). 



References 

Aarseth, S. J. 1999, PASP 111, 1333 

Arnould, M., Goriely, S., & Jorissen, A. 1999, A&A, 347, 572 

Baumgardt, H. & Kroupa, P 2007, MNRAS, 380, 1589 

Baumgardt, H., Kroupa, P, & Parmentier, G. 2008, MNRAS, 384, 1231 

Bedin, L. R., Piotto, G., Anderson, J., et al. 2004, ApJ, 605, L125 

Bekki, K., Campbell, S. W., Lattanzio, J. C, & Noms, J. E. 2007, MNRAS, 267 

Bellini, A., Piotto, G., Bedin, L. R., et al. 2009, A&A, 493, 959 

Bonifacio, P, Pasquini, L., Molaro, P, et al. 2007, A&A, 470, 153 

Brown, J. A. & Wallerstein, G. 1993, AJ, 106, 133 

Caloi, V. & D'Antona, F. 2005, A&A, 435, 987 

Caloi, V. & D'Antona, F 2007, A&A, 463, 949 

CaiTetta, E. 2006, AJ, 131, 1766 

Carretta, E., Bragaglia, A., Gratton, R., & Lucatello, S. 2009, A&A, 505, 139 
Carretta, E., Bragaglia, A., Gratton, R. G., et al. 2006, A&A, 450, 523 
CaiTetta, E., Bragaglia, A., Gratton, R. G., & Lucatello, S. 2008, ArXiv e-prints 
CaiTetta, E., Bragaglia, A., Gratton, R. G., Lucatello, S., & Momany, Y. 2007, 

A&A, 464, 927 
Cottrell, P L. & Da Costa, G. S. 1981, ApJ, 245, L79 

de Mink, S. E., Pols, O. R., Langer, N., & Izzard, R. G. 2009, A&A, 507, LI 
Decressin, T., Baumgardt, H., & Kroupa, P 2008, A&A, in press, 492, 101 
Decressin, T., Charbonnel, C, & Meynet, G. 2007a, A&A, 475, 859 
Decressin, T., Charbonnel, C, Siess, L., et al. 2009, A&A, 505, 727 
Decressin, T., Meynet, G., Charbonnel, C, Prantzos. N., & Ekstrom, S. 2007b, 

A&A, 464, 1029 
Denissenkov, P A. & Herwig, F 2003, ApJ, 590, L99 

D'Ercole, A., Vesperini, E., D'Antona, F, McMillan, S. L. W., & Recchi, S. 

2008, MNRAS, 1228 
Ergma, E. & van den Heuvel, E. P J. 1998, A&A, 331, L29 
Fenner, Y., Campbell, S., Karakas, A. I., Lattanzio, J. C, & Gibson, B. K. 2004, 

MNRAS, 353, 789 

Fraix-Burnet, D., Davoust, E., & Charbonnel, C. 2009, MNRAS, 398, 1706 
Fryer, C. L. 1999, ApJ, 522, 413 

Georgy, C, Meynet, G., Walder, R., Folini, D., & Maeder, A. 2009, A&A, 502, 
611 

Gratton, R., Sneden, C, & Carretta, E. 2004, ARA&A, 42, 385 
Gratton, R. G., Lucatello, S., Bragaglia, A., et al. 2007, A&A, 464, 953 
HaiTis, W. E. 1996, AJ, 112, 1487 

Heger, A., Fryer, C. L., Woosley, S. E., Langer, N., & Hartmann, D. H. 2003, 

ApJ, 591, 288 
Herwig, F 2004a, ApJ, 605, 425 
Herwig, F 2004b, ApJS, 155, 651 

Hirschi, R., Meynet, G., & Maeder, A. 2004, A&A, 425, 649 

Hurley, J. R., Pols, O. R., & Tout, C. A. 2000, MNRAS, 315, 543 

Karakas, A. I. & Lattanzio, J. C. 2003, Publications of the Astronomical Society 

of Australia, 20, 279 
Kobulnicky, H. A. & Skillman, E. D. 1997, ApJ, 489, 636 
Kraft, R. P 1994, PASP 106, 553 
Kroupa, P 2001, MNRAS, 322, 231 
Lada, C. J. & Lada, E. A. 2003, ARA&A, 41, 57 

Lind, K., Primas, F, Charbonnel, C, Grandahl, F, & Asplund, M. 2009, A&A, 

503, 545 
Maeder, A. 1992, A&A, 264, 105 

Marks, M., Kroupa, P, & Baumgardt, H. 2008, MNRAS, 386, 2047 



T. Decressin et al.: Evolution of two stellar populations in globular clusters 



11 



Marks, M., Kroupa, P., & Baumgardt, H. submited, MNRAS 

Milone, A. P., Piotto, G., King, I. R., et al. 2010, ApJ, 709, 1183 

Palous, J., Wunsch, R., Tenorio-Tagle, G., & Silich, S. 2009, in lAU 

Symposium, Vol. 254, lAU Symposium, ed. J. Andersen, J. Bland-Hawthorn, 

& B. Nordstrom, 233-238 
Paresce, F. & De Marchi, G. 2000, ApJ, 534, 870 
Parmentier, G. & Fritze, U. 2009, ApJ, 690, 1 1 12 
Pasquini, L., Bonifacio, P, Molaro, P, et al. 2005, A&A, 441, 549 
Pasquini, L., Bonifacio, P, Randich, S., et al. 2007, A&A, 464, 601 
Pflamm-Altenburg, J. & Kroupa, P 2009, MNRAS, 397, 488 
Piotto, G., Bedin, L. R., Anderson, J., et al. 2007, ApJ, 661, L53 
Piotto, G., Villanova, S., Bedin, L. R., et al. 2005, ApJ, 621, 777 
Portegies Zwart, S. F, Verbunt, F, & Ergma, E. 1997, A&A, 321, 207 
Prantzos, N. & Charbonnel, C. 2006, A&A, 458, 135 
Prantzos, N., Charbonnel, C, & Iliadis, C. 2007, A&A, 470, 179 
Renzini, A. 2008, MNRAS, 391, 354 
Smith, G. H. 2006, PASP 118, 1225 

Tenorio-Tagle, G., Wiinsch, R., Silich, S., & Palous, J. 2007, ApJ, 658, 1 196 

Townsend, R. H. D., Owocki, S. P, & Howarth, I. D. 2004, MNRAS, 350, 189 

Ventura, P & D'Antona, F 2005a, A&A, 431, 279 

Ventura, P & D'Antona, F 2005b, A&A, 431, 279 

Ventura, P & D'Antona, F 2()05c, A&A, 439, 1075 

Ventura, P & D'Antona, F 2()08a, A&A, 479, 805 

Ventura, P & D'Antona, F 2008b, MNRAS, 385, 2034 

Ventura, P & D'Antona, F 2009, A&A, 499, 835 

Ventura, P, D'Antona, F, & Mazzitelli, I. 2002, A&A, 393, 215 

Ventura, P, D'Antona, F, Mazzitelli, I., & Gratton, R. 2001, ApJ, 550, L65 

Villanova, S., Piotto, G., King, I. R., et al. 2007, ApJ, 663, 296 

Wallerstein, G., Leep, E. M., & Oke, J. B. 1987, AJ, 93, 1137 

Wiinsch, R., Tenorio-Tagle, G., Palous, J., & Silich, S. 2008, ApJ, 683, 683 



